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Abstract 

Direct searches for dark matter have prompted in recent years a great deal of excitement within 
the astroparticle physics community, but the compatibility between signal claims and null results of 
different experiments is far from being a settled issue. In this context, we study here the prospects 
for constraining the dark matter parameter space with the next generation of ton-scale detectors. 
Using realistic experimental capabilities for a wide range of targets (including fluorine, sodium, 
argon, germanium, iodine and xenon), the role of target complementarity is analysed in detail 
while including the impact of astrophysical uncertainties in a self-consistent manner. We show 
explicitly that a multi-target signal in future direct detection facilities can determine the sign of 
the ratio of scalar couplings f n /fp, but not its scale. This implies that the scalar-proton cross- 
section is left essentially unconstrained if the assumption f p ~ /„ is relaxed. Instead, we find that 
both the axial-proton cross-section and the ratio of axial couplings a n /a p can be measured with fair 
accuracy if multi-ton instruments using sodium and iodine will eventually come online. Moreover, 
it turns out that future direct detection data can easily discriminate between elastic and inelastic 
scatterings. Finally, we argue that, with weak assumptions regarding the WIMP couplings and 
the astrophysics, only the dark matter mass and the inelastic parameter (i.e. mass splitting) may 
be inferred from the recoil spectra - specifically, we anticipate an accuracy of tens of GeV (tens of 
keV) in the measurement of the dark matter mass (inelastic parameter). 
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I. INTRODUCTION 



Although extensive evidence supports the presence of large quantities of dark matter 
(DM) in the Universe, its specific nature remains undisclosed. Over the last decades, differ- 
ent, complementary techniques have been put forward in the effort to identify dark matter 
and constrain its properties. In particular, direct detection appears very promising in the 
near future and lies today at the cutting edge of research in particle physics and astrophysics. 
The idea behind direct detection is rather simple and compelling: if the dark matter in our 
Galaxy is composed of Weakly Interacting Massive Particles (WIMPs), then the Earth is 
continuously crossing a large flux of WIMPs, which can be detected by their scattering off 
nuclei in underground instruments. The past few years have seen a great excitement in direct 
dark matter searches, mainly due to a mix of both startling signals and null results. On the 
one hand, the DAMA collaboration reports an annual modulation of WIMP-like events with 
impressive statistical significance [Hl2], whereas CoGeNT detects excess events towards the 
threshold energy [3] also with an annual modulation [3]. Under standard assumptions the 
two results are not compatible, but it is worth noting that there is enough freedom on the 
particle physics, astrophysics and experimental sides in order to explain both results with the 
same WIMP candidate (see e.g. [oTflO]). On the other hand, there is a whole array of direct 
detection experiments reporting no events above the expected background. Among them, 
XENON10/100 [TTHT3] and CDMS [HUTU] are perhaps the most sensitive and report small 
numbers of WIMP-like scatterings largely compatible with the background. Also CRESST 
[T7] and EDELWEISS [H] see events in their signal region, but an interpretation of the latter 
in terms of WIMP scattering appears rather difficult. A way out of this puzzling situation 
may be provided by the upcoming generation of direct detection experiments which feature 
ton-scale or multi-ton-scale targets. Hopefully, these instruments will shed light on the topic 
over the next decade and confirm or rule out the dark matter hypothesis. 

Recently, several works have appeared that discuss the prospects for measuring different 
WIMP-nucleon couplings with upcoming direct detection data. For instance, in Ref. [H] (see 
also [201I2T]) the authors address the possibility of extracting the scalar-proton, axial-proton 
and axial-neutron cross-sections with the total event rates in future fluorine, germanium and 
iodine targets. In Refs. [22], [23] (see also [231 [25]) a different method is presented to estimate 
scalar and axial couplings as well as the ratios of proton to neutron couplings. These and 
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other works rely on a certain degree of theoretical assumptions regarding the WIMP-nucleon 
couplings - e.g. one usually assumes that the scalar cross-section dominates over the axial 
cross-section (or vice-versa) or that scalar-proton and scalar-neutron couplings are equal. 
Although these are valid assumptions for the most studied WIMP candidates such as the 
supersymmetric neutralino (see [26]), we should keep in mind that other possibilities exist 
and lead to a rich phenomenology. For example, it is possible to construct models where the 
WIMP scatters off protons and neutrons with different amplitudes - specific models have 
been presented in [HlElEZj- This implies that spin-independent rates no longer scale as A 2 
(A being the target mass number) as usually assumed. Also, as pointed out in the literature 
recently [211 E3 I28ti37] but also ten years ago [3B1 EH], the unknown WIMP local density 
and velocity distribution translate into a sizeable "astrophysical" uncertainty that hinders 
the determination of WIMP properties and has not been addressed in works such as [T9] 1 . 
In this context, the aim of the present paper is to identify what kind of information can 
really be measured in a robust manner with the next generation of ton-scale direct detection 
experiments. We focus our attention on the extraction of WIMP mass, couplings (scalar, 
vector, axial) and inelastic parameter with special emphasis on the importance of target 
complementarity. Ours is an analysis that uses realistic upcoming experimental capabilities 
and keeps particle physics theoretical assumptions to a minimum. We show explicitly the 
effect of relaxing certain widely used assumptions (e.g. the equality of scalar-proton and 
scalar-neutron couplings) on the constrained WIMP parameter space; this issue has been 
ignored in almost all past studies. Moreover, the astrophysical uncertainty is folded self- 
consistently in our results and its role is discussed in detail. We perform a Bayesian analysis 
which allows the translation of statistical and astrophysical uncertainties into the derived 
constraints, thus improving upon analytical approaches as the one presented in [19] . Another 
point addressed here is the possibility of distinguishing elastic and inelastic WIMPs and 
measuring the corresponding inelastic parameter. This study assumes particular relevance 
now since complementarity analyses including searches at the Large Hadron Collider will 
soon be possible, and hence it is essential to know what can - and, perhaps more importantly, 
what cannot - be learnt from direct detection data alone. 

1 Let us stress that other methods [2H [531 123 133 122] have also been developed to extract WIMP properties 
independently from astrophysical unknowns. 
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The paper is organised as follows. In Section [IT] we briefly review the formalism of direct 
searches before specifying the experimental capabilities and dark matter benchmarks in 



Section |III| Section IV outlines our Bayesian methodology while in Section M we present 



the results. Finally, the main conclusions are drawn in Section VI 



II. DIRECT DETECTION FORMALISM 

The central quantity in direct detection studies is the differential rate dR/dE R at which 
WIMP-nucleus scattering events occur in an underground instrument composed of nuclei 
N(A, Z) (for reviews see [4*0l - l4"2"] ). This rate is simply given by the convolution of the WIMP 
flux as seen from the Earth and the scattering cross-section, and it can be conveniently 
written as 

^(Er) = x [*™gF% D {E R ) + ct^FUEr)] X jF{Vmin{E R )) (1) 

in units of counts/ton/yr/keV and where po is the local WIMP mass density, m x is the WIMP 
mass, /i at is the WIMP-nucleus reduced mass and cr fjy (cr i'j^) is the zero- momentum spin- 
dependent (spin- independent) cross-section. The form factors Fsd an d Fgi account for the 
energy dependence of the respective cross-sections. Finally, the factor J 7 is the WIMP mean 
inverse velocity and reads 

F{y min ) = d v , (2) 



v>v. 



rn in. 



f being the local WIMP velocity distribution and v e the Earth velocity in the galactic rest 
frame. The minimum WIMP velocity that produces a nuclear recoil of energy Er is 

where is the nucleus mass and 6 is the inelastic parameter that vanishes in the case of 
elastic scattering and reads 5 = m x > — m x in the case of inelastic dark matter models (x' 
being the excited DM state). 

In equation (JIJ, the astrophysical dependence is encoded in p and J 7 , while the parti- 
cle and nuclear physics enter in the middle factor. Regarding astrophysics, we follow the 
approach of Refs. [35l H3] and model the WIMP velocity distribution as 

f(w) = { \ \ kv o ) ) , ( 4 ) 

for w > v Psr 




where Nf is a suitable normalisation, A; is a shape parameter, v esc is the local escape velocity 
and Vq traces the velocity dispersion. Different parameterisations for the WIMP phase 
space distribution (including dark disks and streams) and the corresponding effect in direct 
detection are studied in j29j [37] . Also, since we are not interested in modelling the annual 
modulation signal, the Earth orbit and the Sun's peculiar motion are disregarded so that v e 
is given by the local circular velocity v^, whose absolute value we shall identify with vq (see 
Ref. [35] for a detailed discussion). 

The WIMP-nucleon scattering cross-section is commonly split into two components, spin- 
dependent (SD) and spin-independent (SI). The former arises from the WIMP-quark axial 
coupling and is expressed as 

^^«^(a p (?) + a^)) 2 , (5) 



in which Gf is the Fermi coupling constant, Jn is the spin of the target nucleus N, a 



p 

(a n ) is the axial WIMP-proton (-neutron) coupling and (Sp) ((S^)) is the expectation 
value of the spin of protons (neutrons) in the nucleus N. The values of (Sp) and (S^) for 
different nuclei may be found in Ref. jH] . Defining the axial WIMP-proton cross-section as 
of D = f /^G|faJ (recall that J p = 1/2, (S£) = 1/2, (S£) = 0), we recast equation @ as 

a x-N = TT — j ° v \\ S p / + \ S n ) a n/a p ) , (6) 

ofip Jn 

where \x p is the WIMP-proton reduced mass. In the body of the paper, we shall use a^ 
and a n /a p as our independent phenomenological parameters. The spin-dependent nuclear 
form factor Fsd is implemented according to the standard prescription [40J: 

FUE R) = f| (7) 

with q 2 = 2ijinE r and 

S(q) = alS 00 (q) + alS n (q) + a Q aiS 01 (q) , (8) 
where a = a p + a n and CL\ — dp — Q"rf 

The structure functions Soo, Sn and Sqi depend on 
the target nuclei and are parameterised following Ref. [15] for the elements used in this work 
(cf . next Section and Table |TJ) . 

Spin-independent scattering arises, instead, from the WIMP-quark scalar and vector 
couplings. The corresponding cross-sections are usually written in terms of the scalar and 



vector WIMP-nucleon couplings f P}Tl and b P)Tl respectively: 



a 



si,o 

X-N 



71 



2 



(Zf p + (A- Z)f n f + — (Zb p + (A- Z)b n ) 



(9) 



Similarly to the SD case, we define the scalar WIMP-proton cross-section a p c = -f-ipfp and 
the vector WIMP-neutron cross-section a^ ec = Tn-Pn^n (the reason for defining a^ ec and not 



a^ ec will become apparent in Sections 



III 



and 



V A ) to recast equation ^ as 



a 



SLO 



x -n ~ -y [< (Z + (A- Z)UIU) + < c (Zb p /b n + (A- Z)f] , (10) 

H'p 



where we made the very good approximation /i n ~ fi p . Again, we will be using the cross- 
sections a p c , <r^ ec and the coupling ratios f n / f p , b p /b n as our phenomenological parameters. 
As for the spin-independent nuclear form factor Fsi, we take the parameterisation introduced 
in [JT] and that has been shown to be a reliable approximation at least for Er < 100 keV 
El: 



where q, r n and s are implicitly expressed in natural units, r 2 = c 2 + ^n 2 a 2 — 5s 2 , s ~ 0.9 
fm, a ~ 0.52 fm and c = (1.23A 1 / 3 - 0.6)fm. 

Now, once one defines a WIMP model (m x , (J p D , a n /d p , cr p c , f n J f p , a^ ec , b p /b n , 6), an 
astrophysical setup (p , v , v esc , k) and a target nucleus N, the differential event rate dR/dE R 
in equation (JI]) is unambiguously fixed. But what direct detection experiments can actually 
observe is the number of events in a given energy bin [E\, E2]: 

/•E2 

N R (E 1 ,E 2 ) = e eff 



dR , . 
dE R —(E R ) 



with 



dR 



dE 



R 



(E 



R) 



dE[ 



R 



2na(E' R ) 



cxp 



dE 



R 



(Er ~ E' R 

2a 2 (E' R 



1 \2 



dR 



dE 



R 



(E\ 



R) 



(12) 



(13) 



where e e // is the mean effective exposure of the experiment and a its energy resolution. The 



number of counts in equation (12) plays a central role in the remaining of this work. 



III. EXPERIMENTAL CAPABILITIES & DARK MATTER BENCHMARKS 



We are interested in assessing the prospects for direct detection experiments over the next 
decade. Hence we focus on a wide range of target elements that are being pushed forward 
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by several collaborations, namely 19 F, 23 Na, 40 Ar, 73 Ge, 127 I and 131 Xe. The respective 
properties and experimental figures are outlined in Table [TTJ In view of the developments 
in the field of direct detection, we consider ton-scale and multi-ton-scale instruments with 
~ O(10) keV threshold energies and energy resolutions parameterised as 

a(E R ) = (a + b^E R /keV + cE R /kevY keV . (14) 



Element 




19p 




eq. (7) of Ref. 05] 


23 Na 




eq. (10) of Ref. 05] 


73 Ge 




eq. (17) of Ref. 05] 


127j 


eq. (20) of Ref. 05 


, "Bonn A" coefficients from Table VI 


131 Xe 


eq. (21) of Ref. 05 


|, "Bonn A" coefficients from Table IX 



TABLE I: The parameterisation of the SD structure functions Soo, Su and Soi used in this work. 



Element 


A Z 


A-Z 




(S%) J N m N [m u ] 


e e// [ton 


.yr] E thr [keV] 


a 


b c d 


19p 


19 9 


10 


0.4751 


-0.0087 1/2 18.998 


1.00 


10 


5.0 


0.0 0.0 1.0 


23 Na 


23 11 


12 


0.2477 


0.0199 3/2 22.990 


1.00 


10 


0.0 


0.245 0.0027 1.0 


40 Ar 


40 18 


22 





39.962 


6.40 


30 


0.0 


0.7 0.0 1.0 


73 Ge 


73 32 


41 


0.030 


0.378 9/2 72.923 


2.16 


10 


0.3 2 


0.0 0.06 2 0.5 


127j 


127 53 


74 


0.309 


0.075 5/2 126.904 


1.00 


10 


0.0 


0.134 0.0008 1.0 


131 Xe 


131 54 


77 


-0.041 


-0.236 3/2 130.905 


2.00 


10 


0.0 


0.6 0.0 1.0 



TABLE II: The nuclear properties of various target elements and the corresponding experimental capabilities. 
The values for (Sp) and (S%) are taken from Ref. [44 and the nuclei mass tum is presented in atomic mass units 
m u = 0.931494 GeV/c 2 . In Section VB we shall use an exposure for 23 Na and 127 I of 3 ton.yr instead of the baseline 



value of 1 ton.yr shown here. See the text for further details. 

Today, two techniques look particularly promising to assemble low-background, ton-scale 
instruments: cryogenic detectors at mK temperatures (using for instance 73 Ge as target ma- 
terial) and noble liquid detectors (featuring for instance 40 Ar or 131 Xe as targets). The former 
technique is being developed and improved by EURECA 07] and SuperCDMS/GEODM 
08], whereas the latter will be used in DARWIN 09], MAX [50] and XMASS [EE], as de- 
tailed in Ref. [35] to which we refer for further discussion. We assume the same experimental 
capabilities for 40 Ar, 73 Ge and 131 Xe as in that paper 2 and summarise the relevant figures 
in Table M 



For the energy resolutions, see in particular Ref. [55] for Xe and Ref. [S3] for Ge. 
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Another interesting target, particularly due its SD sensitivity, is 19 F, already at use in 
experiments as COUPP [Ml 155] . We take an optimistic (but presumably realistic) exposure 
of e e ff( l9 F) = 1 ton.yr and a threshold energy E t h r ( l9 F) = 10 keV. Although COUPP cannot 
in principle yield spectral information, we assume a constant energy resolution of <r( 19 F) = 5 
keV, in line with that obtained in NaF bolometers |56j . 

At last, in order to fully explore target complementarity, also 23 Na and 127 I will be consid- 
ered. These are elements employed in different instruments including DAMA/LIBRA [H[2]; 
iodine is also present in the CF 3 I target material of COUPP [Ml ES]- Similarly to fluorine, 
we consider baseline exposures e e //( 23 Na) = e e jj( 127 I) = 1 ton.yr and threshold energies 
E thr ( 23 Na) = E thr ( 127 l) = 10 keV. For both materials we take a DAMA-like resolution [57], 



a(X = Na,I) = (0.448 ^q x E R /keV + 0.0091 q x E R fkeV) keV (where q x is the quenching 
factor of target X; q^ a = 0.3 and qi = 0.09). 

For all targets presented in Table [IT], we set a maximum recoil energy of 100 keV. It 
is worth noticing that higher energies can lead to tighter WIMP constraints as studied in 
[37] . but it is still not clear if the usual nuclear form factors are valid and if experimental 
background is an issue in that energy range. 

The last ingredient to specify before proceeding with the analysis is the set of DM bench- 



marks to study. We take the WIMP models shown in Table [TTT] that feature various coupling 
configurations. In particular, the SD and SI cross-sections are fixed to values below the 
latest upper limits, 10 -5 pb [5B] and 10 -9 pb [T3] respectively. The reason why we focus on 
50 GeV WIMPs is because most direct detection experiments have an optimal sensitivity 
around this mass - we are indeed interested in extracting the maximum possible information 
from future data and, in that respect, ours is an optimistic work. However, WIMP masses 
of 25 and 250 GeV will also be extensively used. In order to determine the accuracy in 
the measurement of the inelastic parameter 5, we shall also consider a DM benchmark with 
5 = 40 keV in the ballpark of the values explored in the literature [55], and 5 = 100 keV in 



light of recent exclusion limits [60H62] • Table IV shows the expected total event number for 



the various DM benchmarks in Table III and targets in Table O considering a recoil energy 



range between the threshold energy and 100 keV. 
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DM benchmark 


m x [GeV 


up IF" 


(],„ / (la-, 


^[Pb] fn/fp 


< ec [pb] 6 P /6„ 


5 [kcV] 


1 


50 


o 


— 1.0 


i n- 

10 


-9 


1 














o 


la 


25 


o 


—1.0 


10 


-9 


1 














o 


lb 


250 


o 


—1.0 


10 


-9 


1 





u 








o 


2 


50 


o 


— 1.0 


10 


-9 


1 





10 








o 


3 


50 


10~ 5 


— 1.0 


i n- 

10 


-9 


1 














o 


3a 


50 


10 -3 


—1.0 


10 


-9 


1 














o 


3b 


50 


10 -5 


0.0 


10 


-9 


1 





u 








o 


3c 


50 


10~ 5 


+1 


10 


-9 


1 














o 


3d 


50 


10~ 3 


0.0 


10 


-9 


1 





u 








o 


3e 


50 


10 -3 


+1 


10 


-9 


1 














o 


3f 


25 


10 -5 


—1.0 


10 


-9 


1 





u 








o 




250 


10~ 5 


—1.0 


10 


-9 


1 














o 


3h 


25 


10 -3 


—1.0 


i n- 

10 


-9 


1 














o 


3i 


250 


10 -3 


—1.0 


10 


-9 


1 














o 


4 


50 


o 


—1.0 


10 


-9 


1 





u 








40 


4a 


50 


o 


—1.0 


1U 


-5 


1 
1 


U 


n 
U 


n 
U 


U 


100 


4b 


25 





-1.0 


io- 


-9 


1 














40 


4c 


250 





-1.0 


10" 


-9 


1 














40 


5 


50 


10~ 5 


-1.0 


io- 


-9 


1 














40 


5a 


25 


10~ 5 


-1.0 


10" 


-9 


1 














40 


5b 


250 


10~ 5 


-1.0 


10" 


-9 


1 














40 



TABLE III: The properties of the DM benchmarks used to assess the prospects of the upcoming generation of 
ton-scale direct detection experiments. 

IV. METHODOLOGY 



Given the formalism presented in Section [XT] and the DM benchmarks and experimental 



capabilities specified in Section III , we start by generating mock ( "true" ) data for each bench- 
mark and each target nucleus assuming our fiducial astrophysical model: po — 0.4 GeV/cm 3 , 
Vq = 230 km/s, v esc = 544 km/s, k — 1. Ten linearly-spaced energy bins are constructed 
between the relevant threshold energy and 100 keV. For the different DM benchmarks and 
target materials, bin counts range from a few to several hundreds. Along the work, we shall 
use the mock data corresponding to a specific target, or a combination of data sets as shown 
in Table |Vj Data I employ nuclei widely spread across atomic number Z; data II feature the 
targets used in Ref. [19]; data III include all targets and, finally, data IV are a variant of 
data III with enhanced 23 Na and 127 I exposures (3 ton.yr). The usefulness of each data set 
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DM benchmark 


19p 


23 Na 40 Ar 


73 Ge 


127j 


131 Xe 


1 


11 


1 o 

18 


i or; 

Izo 


303 


183 


362 


la 


10 


lb 


o o 
Z5 


161 


71 


142 


lb 


3 


b 


Q 

oo 


145 


70 


159 


2 


14 


zl 


1 CO 

163 


399 


245 


488 


3 


282 




i or; 

Izo 


387 


193 


378 




26975 


1 1 f\ A A 

111)44 


i or; 

125 


8761 


3548 


1955 


3b(*) 


271 


1 on 

180 


i or; 

125 


305 


603 


371 


3cf*l 


262 


203 


i or; 

lz5 


434 


634 


406 


3df*l 


26114 


13046 


i or; 

lz5 


454 


5863 


1053 


3ef*) 


25168 


15zlb 


i or; 

125 


13330 


8887 


4549 


3fY*) 


284 


1/11 

141 


28 


211 


224 


148 


3z(*) 
°t>\ ) 


101 


1 

bl 


o o 
Ob 


176 


260 


168 




27303 


ybyi 


OQ 

2b 


5113 


1374 


758 




9594 


4251 


o o 

88 


3626 


1731 


977 


4 


<1 


<i 


i 

4 


28 


31 


64 


4a 


<1 


<i 


<1 


<1 


425 


1025 


4b 


<1 


<i 


<1 


2 


2 


4 


4c 


<1 


<i 


14 


40 


38 


72 


5 


<1 


<i 


4 


38 


33 


66 


5a 


<1 


<i 


<1 


3 


2 


4 


5b 


<1 


<i 


14 


52 


37 


75 



TABLE IV: The expected total number of events in each of the targets in study and for the various DM benchmarks. 
The recoil energy range considered spans from the threshold energy up to 100 keV. The experimental capabilities 
assumed are the ones in Table |TT| except for the lines marked with (*) in which an exposure of 3 ton.yr was taken for 
23 Na and 127 I. Notice that even in the cases where a very small number of events (<1) is expected, the corresponding 
target data is valuable in constraining the WIMP parameter space. 

will become apparent in the next Section. 



Data set 


targets comments 


data I 


19 F+ 40 Ar+ 73 Ge+ 131 Xe 


data II 


19 F+ 73 Ge+ 127 T 


data III 


19 F+ 23 Na+ 40 Ar+ 73 Ge+ 127 I+ 131 Xe _ 


data IV 


19 F+ 23 Na+ 40 Ar+ 73 Ge+ 127 I+ 131 Xe £e// (23 Na) = ^(127^ = 3 ton . yr 


TABLE V: The mock experimental data sets used along the work. 



The next step is to scan over the parameter space composed by the WIMP properties 
(m x , <7p D , a n /a p , a s p c , f n /f p , < ec , b p /b n , 5) and the astrophysical unknowns (p Q , v , v esc , 
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k), whose ranges are outlined in Table VI To this end, we use the MultiNest code [63H65] 3 
which is a very efficient sampler of higher-dimensionality parameter spaces and relies on 
Bayes' theorem to update the prior p(9): 

where p(9\d) is the posterior probability function, 9 (d) represents the parameter (data) 
set, C{9) is the likelihood of the parameter set 9 given the data d and p(d) is the so-called 
Bayesian evidence (for a review on Bayesian tools and techniques see Ref . [66] ) . In this work 
the MultiNest code is run with nlive=3000 live points, an efficiency parameter |ef f =0"T3 and 



a tolerance |tol=0 .8| - these parameters were found to be appropriate to scan efficiently and 
in a reasonable period of time the multi-dimensional parameter space in study; see [SZ] for 
further details. Since we are interested in parameter estimation, the Bayesian evidence is 
merely a normalisation constant and will be dropped in the following. We take the priors 



specified in Table VI and implement a binned likelihood defined as 



^)=nnTrT ex p(- iv ^) > 

% b ' 

in which i indexes the experimental setup associated to each target material, b indexes 
the (ten) energy bins in each experimental setup, is the number of counts in bin b and 



target i for the parameter set 9 (computed using equation ( 12 )) and is the corresponding 
number in the mock data (also computed with equation ( 12 ) but taking the true parameter 
set). Several comments are in order here. Firstly, we assume negligible background events, 
which is optimistic but appears reasonable in light of the prospects for low-background, 
ton-scale direct detection experiments [3ZJ HHl EE]. It would be interesting to study the 
impact of background in the findings presented here and, in particular, to understand how 
target complementarity can help interpreting a WIMP signal in the presence of background. 
A detailed study including background modelling is left for future work. Secondly, we 
use the binned likelihood because it is perhaps the most straightforward way to derive 
constraints without the need to worry about realisation noise. In other words, our "true" 
counts are directly drawn from the corresponding event rate without including Poisson 



3 The author acknowledges Louis Strigari and Roberto Trotta for kindly providing access to their private 
direct detection code. 
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scatter. Therefore our results do not include possible realisation noise; this effect has been 
included in other works in the literature [331 [Ml EZ] where the unbinned likelihood was 
implemented. 



1—* O Y* O TT1 /"\T" Ay 

JTcxI dint; IfcJI 


i aiigt; 


pi lOI 




nCLUCldl 


log 10 (m x /( J e\ ) 


1 Q nl 

o.U J 


lldl 




lnrr ^0 lncr 9^ lnrr 9^0 
1O §10 0U ' iO 610 /0 ' iO &10 Z0U 
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TABLE VI: The phenomenological WIMP-related parameters considered to explore the capabilities of the next 
generation of direct detection instruments. The second and third columns specify the range and prior fed into the 
MultiNest code, while the last column indicates the fiducial values. 



Regarding the astrophysical model - encoded in the parameters p , v , v esc and k three 
setups will be used: (i) "fixed" astrophysics, where we set the astrophysical parameters to 



their fiducial values shown in Table VI (ii) "varying" astrophysics, where uncertainties on 
the local DM density and velocity distribution are taken into account following Ref. [35J, 
namely a flat prior on the shape parameter k and la priors on po, ^o and v esc according 



to Table VI and (Hi) "flat" astrophysics, where po ; v esc and k are all varied in their 
ranges with flat priors. The first approach has been extensively used in the literature and 
assumes an exact knowledge of the astrophysical setup, which is certainly not the case as 
of today. The second setup, instead, is a reasonable assessment of present uncertainties 
on the different astrophysical unknowns (see |35j for a detailed discussion), while the third 
approach assumes rather poor astrophysical knowledge and is useful to study the constraints 
that can be placed on the parameters p$, Vq, v esc and k with direct detection alone. 

Let us finally point out that, since our phenomenological parameters include cross-sections 
(o~p D , <7p C , cx^ ec ) and WIMP-nucleon couplings (a n /a p , f n / f p , b p /b n ) - instead of WIMP-quark 
couplings -, we do not fold in nuclear uncertainties pertaining the spin content of nucleons 
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and the 7r-nucleon sigma term. For alternative analyses done in the framework of specific 
WIMP models, see [2S1E3]. 



V. RESULTS 

A. Spin-independent couplings 

The first question we wish to address is whether ton-scale direct detection experiments will 
be able to determine scalar-proton and scalar-neutron couplings independently, i.e. measure 
both <7p C and /„/ f p in equation ( JlO| . Notice that although the widely studied supersymmetric 
neutralinos feature f p ~ f n [25], it is rather easy to build WIMP models with uncorrelated 
proton and neutron couplings f p and f n - for specific implementations see [6j d |27] . This 
means that the WIMP-proton and WIMP-neutron scattering amplitudes may interfere non- 
coherently and the well-known scaling dR/dEft oc A 2 for spin-independent searches is no 
longer valid. For the time being vector and axial couplings are neglected. At fixed DM 
mass, the normalisation of the recoil spectrum is basically set by the quantity 

I sc = a s p c (Z + (A-Z)f n /f p ) 2 . (17) 

Therefore, single-target data cannot break the degeneracy between a p c and f n J f p . Instead, 
by using targets with distinct ratios (A — Z)/Z, one can hope to break the degeneracy and 
constrain scalar-proton and scalar-neutron couplings. However, for heavy nuclei - the most 
sensitive to SI scattering - nuclear stability entails an almost universal ratio (A — Z)/Z ~ 1.4. 
That is why the authors of Ref. [25] suggested the use of target elements with Z < 17 (for 
which (A — Z)/Z ~ 1) along with heavy materials in order to measure f n /fp- 



Here, we test such suggestion by taking the DM benchmark 1 in Table HI as the true 
model and focussing on 19 F, 40 Ar, 73 Ge and 131 Xe as target materials. Notice that the targets 
are chosen across the periodic table to maximise their complementarity regarding the ratio 



(A — Z)/Z. Following the procedure outlined in Section IV, we use the mock data to set 
constraints on the parameter space defined by m x , a p c and f n / f p (all other parameters are 
kept fixed to their input values corresponding to benchmark 1 in Table III). Figure [l] shows 



the joint 95% posterior probability contours derived with 40 Ar (left frame, grey), 131 Xe (left 
frame, yellow) and 19 F+ 40 Ar+ 73 Ge+ 131 Xe (right frame, red) keeping the astrophysical pa- 



rameters fixed. The solid (dashed) line in the left frame is an iso-contour of the quantity (17) 
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FIG. 1: The joint 95% posterior probability contours for the DM benchmark 1 and fixing the astrophysical pa- 
rameters. For simplicity the joint 68% posterior probability contours are not shown. In the left frame, the grey 
and yellow contours show the reconstruction capabilities of 40 Ar only and 131 Xe only respectively, while in the right 
frame the red contour corresponds to 19 F+ 40 Ar+ 73 Ge+ 131 Xe. The solid (dashed) curve in the left frame represents 
the degeneracy direction according to equation ( 171 for 40 Ar ( 131 Xe). As the Figure shows, a multi-target signal can 
determine the sign of f n /f P , but not its absolute value. 

in the case of a 40 Ar ( 131 Xe) target. It is evident from Figure [I] that target complementarity 
is particularly effective for determining the correct sign of the /„/ f p ratio. Nevertheless, dif- 
ferent targets constrain essentially the same parameter space for positive values of f n /fp, arid 
the degeneracy between a p c and f n j f p cannot be entirely lifted. This degeneracy leads to a 2a 
uncertainty on a p c of about 1.5 orders of magnitude (i.e. 8 x 10~ n pb < a p c <6x 1CT 9 pb), 
as can be better appreciated from Figure [2] where the effect of marginalising over the as- 
trophysical parameters is shown by the green contours. Notice that the flat bottom-end of 
the contours in the right plot of Figure [2] would extend to lower o s v c if the prior on f n /f p 



was not restricted to values smaller than 4. The prior on the ratio f n /fp (check Table VI) 



was intentionally stretched beyond the values found in the literature in order to assess what 
direct detection experiments can really tell us about this parameter independently of the- 
oretically motivated assumptions. An even more conservative approach would be to adopt 
a log prior on f n / f p (that assumes no knowledge about the scale), but we feel that the flat 
prior on this parameter is enough for the purposes of the present work. Let us stress at this 
point that the poor accuracy on the determination of a s p stems precisely from the relaxation 
of the hypothesis f n /fp = 1 (and not from astrophysical uncertainties). Of course, if one 
is interested in DM candidates with ~ f p such as supersymmetric neutralinos, then the 
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FIG. 2: The joint 95% posterior probability contours for the DM benchmark 1 and the data set 
19 F+ 40 Ar+ 73 Ge+ 131 Xe. The red contours are the same as in Figure [l] (right), while the green contours show 
the effect of marginalising over the astrophysical uncertainties. The vertical dashed lines in the left frame indicate 
the range 0.8 < fn/fp % 1.2. Without a priori assumptions on the ratio fn/fp, the next generation of direct detec- 
tion instruments can only pinpoint o s ^ within 1.5 orders of magnitude (if m x ~ 50 GeV). The right plot also shows 
the 95% posterior probability contours for the DM benchmarks la and lb when marginalising over astrophysical 
uncertainties. 

contours are restricted to 0.8 < fn/fp ^ 1-2 (including nuclear uncertainties, see [26j ) , as 
shown by the vertical dashed lines in Figure [2] (left), and the accuracy on a s v c is drastically 
improved to about 0.5 orders of magnitude (i.e. 5 x 10~ 10 pb < a s p c < 2 x 10 -9 pb) when 
marginalising over the astrophysical parameters. We have checked that the contour tail at 
negative values of f n j f p reduces drastically (but does not entirely disappear) when using 
the target combination 19 F+ 23 Na+ 40 Ar+ 73 Ge+ 127 I+ 131 Xe (data III) and marginalising over 
the astrophysical unknowns. 

On the other hand, an interesting result of our analysis is the fact that the DM mass 
can be tightly constrained (if m x ~ 50 GeV) despite of the a s p c — /„/ f p degeneracy and the 
accuracy attained is in line with the ones found in works assuming f n = f p (see e.g. [55] 1ST]). 
Notice that this statement does depend on the true WIMP mass: for comparison, we show 
on the right plot of Figure [2] the constraints for the case of m x = 25 GeV (benchmark la) 
and m x = 250 GeV (benchmark lb). The mass accuracy for the three fiducial masses is 
similar (slightly better, in fact) to the findings of Ref. [35], even though we use here one 
extra target ( 19 F) and let f n /f p vary. In particular, the 95% lower limit on the DM mass 
for the case m x = 250 GeV shown in Figure p] (right) is slightly more stringent than in 
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Ref. [35] given the use of 19 F here that was not considered in that work. The accuracy on 
<jp C , instead, is much worse than in Ref. J3H] and that is due to the a^ c — f n / f p degeneracy 
explained above. Finally, let us notice that astrophysical uncertainties - whose impact is 
stressed in Figure [2] - hinder the determination of all parameters, but leave qualitatively 
unchanged all the points discussed above. 

Overall, the punchline of this exercise is two-fold. Firstly, the next generation of ton- 
scale direct detection experiments will be able to determine the sign - but not the scale - 
of fn/ fp- This means it will be virtually impossible to tell apart DM candidates with pure 
scalar-proton and pure scalar-neutron couplings. However, the upcoming instruments have 
the capability to test isospin-violating dark matter [6j [7], [91 QUI E7] as an explanation of 
DAM A/LIBRA [H E] and CoGeNT [3] observations compatible with XENON10/100 [T[> 
[13] and CDMS [HI EE] upper limits - recall that in Ref. [Zj a value f n J f p ~ —0.7 was used 
and this is falsifiable in light of Figure [TJ Let us note that the constraints presented in 
Figures [l] and [2] may vary significantly depending on the true value of /„/ f p ; a full analysis 
using several true values for this coupling ratio is beyond the scope of this paper but is 
an interesting topic for future work. Secondly, it seems evident from Figure [2] (right) that 
direct detection alone is able to extract the DM mass precisely (if m x ~ 50 GeV) but shall 
not determine a s v c to better than one order of magnitude. This is an important guideline for 
future complementarity studies including both direct detection and accelerator searches. 

Up to now we have assumed vanishing vector couplings. An interesting question is to 
know whether scalar-interacting and vector-interacting DM candidates can be distinguished 
with future direct detection data. Since Majorana particles do not have vector couplings 
but Dirac particles do, this would be a first step in the determination of the nature of dark 



matter. Clearly, according to equation (10), a DM candidate featuring scalar couplings only 
ip n = bp = 0) with a given f n / f p produces exactly the same recoil spectrum as a particle of 
the same mass and cross-section with vector couplings only (f n = f = 0) and b' n /b' = f n j ' f p . 
These scenarios cannot be told apart. Instead, we shall focus on candidates with f n /fp = 1 
and b p /b n = 0. The former is a very typical value for supersymmetric neutralinos (see [26]). 
while the latter is featured by sneutrinos and heavy neutrinos [10]. In this case, target 
complementarity may play a role in extracting both scalar and vector couplings since the 
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FIG. 3: The joint 95% posterior probability contours for the DM benchmark 2 and the data set 
i9 F+ 40 Ar+ 73 Ge+ i3i Xe Thg red contours 

were obtained for fixed astrophysics, while the green contours include 
astrophysical uncertainties. Even disregarding astrophysical uncertainties, the distinction between scalar and vector 
scatterings is virtually impossible. 

rate normalisation (for a given mass m x ) scales as 

Lc+vec = a s p c A 2 + a v :\A - Zf (f n /f p = l,b p /b n = 0) . (18) 

Now, adopting the DM benchmark 2 in Table |III| and following the same procedure as 
before, we scan over the parameter space (m x , at 10 , cr^ ec ) while fixing f n / f p = 1, b p /b n = and 
all other parameters to their true values. The joint 95% posterior probability contours are 
shown in Figure[3]for 19 F+ 40 Ar+ 73 Ge+ 131 Xe (data I) and the cases of fixed (red) and varying 
(green) astrophysics. As evident from this plot, the prospects for distinguishing between 
pure scalar, pure vector or mixed DM particles are rather poor. Notice, in particular, that 
no constraint can be placed on a p c nor a^ ec within their priors - and this statement holds 
even without marginalising over astrophysical uncertainties. As in the scalar-only scenario, 
different true values for f n / f p and b p /b n may lead to different results but such study is 
deferred to future work. 

B. Spin-independent and spin-dependent couplings 

We now turn to the prospects of measuring SI and SD cross-sections with upcoming direct 
detection experiments. Despite the a p c — f n / ' f p degeneracy discussed in the previous Section, 
here we shall fix f n j f p = 1 and work in the parameter space (m x ,a p c ,a p D ,a n /a p ). Once 
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again, vector couplings are neglected. This has been the setup used in several works such 
as [T5] or [33]. As pointed out in Ref. [19] and clear from equations ^ and (10) in Section 



PH the different couplings can be constrained by using a target with good SI sensitivity and 
targets with distinct (S p )/(S^) in order to lift the degeneracy between and a n /a p . 
Indeed, while the normalisation of the SI rate in a given material (for fixed m x ) is given by 



equation (17), the SD rate depends on 

Isd = * S P D {{S?) + (S»)a n /a p ) 2 . (19) 



The iso-contours of this quantity for different nuclei are shown in Figure [4] (top left and 
bottom left frames) in the case of benchmark 3. The complementarity between the targets 
is clear and that is the key to measure simultaneously <jp D and a n /a p as we shall see. Notice 
that the DM benchmark 3 features a n /a p = —1 which is a typical value when considering 
supersymmetric neutralinos - in fact, these usually present a n > and a p < [191 133] . 

To begin with, we consider the targets proposed in [15] . i.e. 19 F+ 73 Ge+ 127 I. Using the 
mock counts associated to benchmark 3 and each target material, we shown in the top plots 
of Figure [4] the joint 95% posterior probability contours for 19 F+ 73 Ge+ 127 I (data II) and 
the cases of fixed (red) and varying (green) astrophysics. Note that the contours are well 
inside the a n /a p range, which justifies our flat prior on this parameter. Now, the posterior 
is markedly bimodal: there is a ("correct") solution for a n /a p < and a ("wrong") solution 
for a n /a p > 0. This is precisely what was found in [19] where, using total event rates, 
the authors identified one solution for each sign of a n /a p . Our analysis, however, makes 
use of the full energy spectrum and includes the impact of marginalising over astrophysical 
uncertainties. Moreover, for fixed astrophysics we anticipate a 2a measurement a n /a p = 
0.85 ± 0.36 (restricted to a n /a p > 0) and a n /a p = —0.95 ± 0.38 (restricted to a n /a p < 0) 
with 19 F+ 73 Ge+ 127 I. Marginalising over the astrophysical parameters results in very similar 
a n /a p constraints as can be seen by the contours in the top left frame of Figure |ij In 
particular, the absolute value of a n /a p can be well constrained, but not its sign. This is 
essentially opposite to the case of f n /fp for which the sign can be determined but not its 
scale, cf. Section V A and Figures [I] and [2] The reason for such difference lies in the different 



signs and values of (A — Z)/Z (for f n / f p ) and (S^)/(Sp) (for a n /a p ) featured by the nuclei 
at use. Indeed, the values of (A — Z)/Z just slightly vary between 1.11 ( 19 F) and 1.43 
( 131 Xe), while the ratio (S^)/ (S p ) may be either positive or negative and either very small 
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FIG. 4: The joint 95% posterior probability contours for the DM benchmark 3. In the top frames we show the 
reconstruction capabilities of 19 F+ 73 Ge+ 127 I for fixed (red) and varying (green) astrophysics. In the bottom frames 
we have marginalised over astrophysical uncertainties and show the effect of using different data sets. The curves 



in the top left and bottom left frames represent the degeneracy directions according to equation ( 19 1 for different 



target nuclei. The "wrong" solution at a„/a p > can be discarded by using all target data with an exposure of 3 
ton.yr for 23 Na and 127 I. 



(—0.018 for 19 F) or very large (12.6 for 73 Ge). Also, notice that a true ratio a n /a p ^ — 1 
leads in principle to significantly distinct constraints on the WIMP parameter space, as we 
shall see below. In any case, despite the degeneracy between a^ and a n /a p , we find that 
m x , <7p C and <r'p D are determined to a good accuracy using simply 19 F, 73 Ge and 127 I for a 
DM particle with the properties of benchmark 3. 

For comparison, we focus on Ref. [33J where a rather complete study in the context of 
supersymmetry was presented. The authors have used the projected capabilities of ton-scale 
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xenon, germanium and COUPP-like experiments to constrain the supersymmetric parameter 
space, including both astrophysical and nuclear uncertainties. For their DM benchmark 1 
(m x ~ 108 GeV, a s p c ~ 3.9 x 10~ 9 pb, a s p D ~ 2.8 x 10~ 5 pb, a n /a p ~ -0.87), they are able 
to measure all three cross-sections (a p c , a p , cr^) and find that the 2a uncertainty on a p 
and a n is of order 50% - check their top right plot in Figure 11. It is reassuring that this 
uncertainty is in line with the values stated in the last paragraph, although the analyses are 
not strictly comparable. Let us notice, however, that the results in Ref. [33] are basically 
restricted to a p < 0, a n > (or, in our notation, a n /a p < 0) because the authors are working 
on the supersymmetric framework, which means that a possible solution at a n /a p > (our 
right blob in the top left plot of Figure [4]) is dismissed a priori. 

In order to test if further direct detection data are able to break the a p D — a n /a p degen- 
eracy in Figure [4] (top frames), we have rerun the scan but now making use of all targets, 
19 F+ 23 Na+ 40 Ar+ 73 Ge+ 127 I+ 131 Xe (data III). The corresponding posterior probability con- 
tours are displayed in grey in the bottom plots of Figure |4| where we have marginalised 
over astrophysical uncertainties. The wrong solution - at a n /a p > - is now more tightly 
constrained, but cannot be excluded at the 2a level. This persistence of the dual solution is 
due to the fact that three target nuclei ( 19 F, 73 Ge and 131 Xe) are approximately degenerate 
in the regions around a n /a p ~ ±1 as shown by the line contours in Figure [4] (bottom left 
frame). The dual-shaped posterior does depend on the true ratio a n /a p (here a n /a p = — 1) 
and its exact behaviour is in principle different for distinct true a n /a p values - that is an 
interesting topic that we leave for future research (cf. also Figure [5]). Eventually, to com- 
pletely discard the posterior volume at a n /a p > 0, one can bet on a larger exposure for 
targets like 23 Na and 127 I that effectively break through the degeneracy. That is what we 
have done by considering a fourth data set including all targets as in data III but assum- 
ing enhanced 23 Na and 127 I exposures, e e //( 23 Na) = e e //( 127 I) = 3 ton.yr. The results are 
shown in the bottom plot of Figure [4] by the cyan contours and entail the 2a measurements 
!ogio(^p D /P b ) = -5-0 ± 0.25 or a^ D = 1.0+g;J| x 10~ 5 pb and a n /a p = -0.97 ± 0.40 when 
marginalising over the astrophysical parameters. In this case, the positive a n /a p solution is 
entirely ruled out (at ~ 10a). This means that both a p D and a^ (or equivalently a p D and 
a n /a p ) are measurable with future direct detection data. 

Up to now we have focussed on the DM benchmark 3 featuring m x = 50 GeV, a p = 10~ 5 
pb and a n /a p = — 1. It is however important to check how our results depend on the true 
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FIG. 5: The joint 95% posterior probability contours for the DM benchmarks 3 and 3a-i using data IV and 
marginalising over the astrophysical parameters. The plots show how the accuracy in the extraction of WIMP 
parameters can change dramatically according to the true underlying WIMP properties. 



DM mass m x and SD coupling ratio a n /a p . Also, the latest experimental limits [58| |69| ITU] 
still allow SD cross-sections as large as 10~ 3 pb. Therefore, the additional DM benchmarks 



3a-i in Table III are studied and the corresponding 95% posterior probability contours using 



data IV and marginalising over astrophysical unknowns are shown in Figure [5j As can be 
seen in the left frame of this Figure, the a n /a p uncertainty is a complicated function of the 
underlying true value of a n /a p and appears larger for true values a n /a p ~ than \a n /a p \ ~ 1 



SD 



10" 



(when a p D = 1CT 5 pb). Moreover, for the DM benchmarks 3a, 3d and 3e with a, 
pb, the uncertainties shrink considerably given that the fraction of SD events relative to SI 
events is much larger than for <jp D = 1CT 5 pb. In the right frame of the same Figure one can 
appreciate the reconstruction capabilities for different values of m x and a p - as expected, 
high masses and low cross-sections lead to larger uncertainties. 



C. Inelastic parameter 

In direct dark matter searches, the WIMP-nucleus scattering is widely assumed to be 
elastic when deriving constraints on the WIMP properties. This basically amounts to setting 
5 = in equation (|3]). However, as long noticed in the literature [59], inelastic dark matter 
scenarios are perfectly viable and lead to a very diversified phenomenology. In fact, while in 
the elastic case (5 = 0) v min oc a/Er and thus the recoil rate rapidly decays with increasing 
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FIG. 6: The joint 68% and 95% posterior probability contours for the data set 19 F+ 40 Ar+ 73 Ge+ 131 Xe and marginal- 
ising over the astrophysical uncertainties. The different contours were obtained assuming the DM benchmark 1, la, 
lb, 4, 4a, 4b and 4c. In the left frame the constraints are adjacent to the corresponding DM benchmark and in the 
right frame the same colour coding is used. The left plot shows that the next generation of instruments will be able 
to determine the inelastic parameter 5 with an accuracy of tens of keV. 



energy, inelastic dark matter (5 > 0) induces a behaviour v m i n ex 1/ a/E/? (v m in ^ VEr) 
at sufficiently small (large) recoil energies which translates into a peak in the event rate 
(damped by the nuclear form factor). Inelastic DM models featuring 5 ~ — 100 keV have 
been proposed |59j to reconcile DAMA/LIBRA findings with the null results from other 
direct detection experiments, although these models are now in tension with present data 
[60H62] . It is therefore an important matter to study the prospects of measuring - or at 
least constraining - the inelastic parameter 8 with future ton-scale instruments. Ref. [32] has 
studied how 8 (and m x ) can be constrained independently of astrophysics. Here, we pursue 
a slightly more aggressive approach by modelling astrophysical uncertainties as outlined in 



Section |IV[ and present the expected accuracy on 5 once the next generation of instruments 
comes online. For definiteness, we assume the DM benchmarks 1, la, lb (5 = 0), 4, 4b, 4c 
(5 = 40 keV) and 4a (5 = 100 keV), while scanning over (m x ,o"p C , 5) and fixing all other 
parameters to their true values. Notice that for the DM benchmark 4a with 5 = 100 keV 
we use an enhanced scalar cross-section a sc = 10~ 5 pb (still perfectly compatible with null 
searches for S = 100 keV and m x = 50 GeV, cf. Figure 4 in Ref. [60J) in order to have a non- 
negligible event rate for the used targets. Figure [6] shows the joint 68% and 95% posterior 
probability contours for the different DM benchmarks when using 19 F-r- 40 Ar+ 73 Ge+ 131 Xe 
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(data I) and marginalising over the astrophysical parameters. It is remarkable that the 
complementarity between different targets enables the measurement of 5 rather accurately. 
In particular, we obtain a 2a range 5 = 4 ± 6 keV for benchmark 1, 5 = 36 ± 17 keV 
for benchmark 4 and 5 = 98 ±11 keV for benchmark 4a, all featuring m x = 50 GeV. 
Moreover, similar results are obtained for the DM benchmarks with m x = 25 GeV (la, 4b) 
and m x = 250 GeV (lb, 4c) as evident from Figure [6j This is a very interesting result 
indeed: even in the case of benchmarks featuring m x = 250 GeV, for which the DM mass 
itself can only be bounded from below, the inelastic parameter 5 can be robustly pinpointed. 
Recall that the results in this Section are obtained assuming f n /f P = 1, so care must be 
taken in interpreting these figures in light of the discussion in Section |V A Finally, we have 



checked that there is a significant correlation between 5 and vq, so that a vo measurement 



tighter than presented in Table [VI] may lead in the future to a better determination of the 
inelastic parameter. 

D. General case 

In order to make the analysis tractable, we have made thus far some simplifying assump- 



tions, namely or = erf = in Section VA and f n / f p = 1 in Sections VB and VC Here, 



we let all parameters in Table VI (except for a^ ec and b p /b n ) vary and identify which are 
the ones that can be robustly pinpointed. Adopting benchmark 5 as the true model and 
using 19 F+ 23 Na+ 40 Ar+ 73 Ge+ 127 I+ 131 Xe (data III), we scan over the entire parameter space 
(check Table |VI[ ) and marginalise over astrophysical uncertainties obtaining the results dis- 
played in Figure [7} It turns out that only two WIMP properties can be measured with good 
accuracy: m x and 5. The respective 2a ranges read m x = 52 ±20 GeV and 5 = 35 ± 15 keV. 
These figures are very robust against astrophysical uncertainties and rely on a minimum 
of theoretical assumptions. Apart from m x and 5, all other parameters are left essentially 
unconstrained within their priors. For completeness we show in the top right and bottom 
plots of Figure [7] the cases of benchmarks 5a (m x = 25 GeV) and 5b (m x = 250 GeV). While 
for 5a both m x and 5 can be reasonably well pinpointed, for 5b only the mass splitting S 
can be measured. It is interesting to notice that for m x = 250 GeV the SD cross-section 
a s p D is systematically underestimated. 

Finally, we repeat the scan for benchmark 5 but this time with a "flat" astrophysical 
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FIG. 7: The joint 95% posterior probability contours for the DM benchmarks 5 (green), 5a (red) and 5b (yellow) 
using data III and marginalising over the astrophysical uncertainties. For clarity the contours corresponding to 
benchmarks 5a and 5b were omitted in the top left frame. The only two parameters that can be well constrained 
within their priors are the DM mass m x and the inelastic parameter <5. 



setup (as defined in Section IV). This is an useful procedure to check the dependence of 
WIMP signals on astrophysical parameters. We find that m x and 5 can still be robustly 
determined, although with larger uncertainties than stated in the previous paragraph. Unfor- 
tunately, however, in this case a multi-target detection with the next generation of ton-scale 
instruments is not enough to constrain the astrophysical parameters po, Vq, v esc and k - in 
fact, the posterior distribution on these quantities is far wider than the priors used for the 



astrophysical uncertainty in Table VI 
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VI. CONCLUSION 



In this work we have fully explored the ability of the next generation of ton-scale direct 
detection experiments to constrain WIMP properties. Using realistic upcoming experimental 
capabilities and including astrophysical uncertainties, we have studied how target comple- 
mentarity can be used to improve WIMP-related constraints in case of a multiple-target 



analysis. Our main findings - outlined in Table VII - may be summarised as follows: 



It will be possible to determine the sign of the ratio of scalar couplings /„/ f p , but not 
its absolute value. This entails an uncertainty on the scalar-proton cross-section o s v c 
of more than one order of magnitude (8 x 10 -11 pb < o s v c <6x 10~ 9 pb), if we are to 
relax the widely used (but seldom tested) assumption f n / f p = 1. 

Scalar and vector cross-sections cannot be isolated even using multi-target data. Thus, 
the prospects for discriminating Majorana and Dirac WIMPs with upcoming ton-scale 
instruments are rather pessimistic. 

The absolute value of a n /a p can be fairly well constrained. However, we found that 
even with ton-scale direct detection experiments, two disconnected regions in the pa- 
rameter space (a n /a p ,<Tp D ) are compatible with the mock data: one at a n /a p < (the 
"correct" solution) and one at a n /a p > (the "wrong" solution). The latter can be 
entirely discarded making use of targets such as 23 Na and 127 I which are particularly 
complementary to 19 F, 73 Ge and 131 Xe. But we could only isolate the solution at 
c^n/cLp < with exposures for 23 Na and 127 I of 3 ton.yr that are probably out of reach 
within the next decade. 

The inelastic parameter 5 can be measured to an accuracy of tens of keV. This state- 
ment holds even for DM masses as large as 250 GeV in which case the mass itself 
cannot be precisely measured. Furthermore, inelastic DM candidates with modest 
mass splittings of S ~ 40 keV can be easily distinguished from the standard elastic 
scattering scenario. 

The only WIMP parameters that can be extracted from the data in a very robust 
manner are the DM mass m x and the inelastic parameter 5. All other parameters can 
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Section data set benchmark 


m v [GeVl 


f If 

J TL / J p 


crl c [10" 9 pbl 

p L 1 J 


a In 


a^ D [10" 5 pb] 5 [keV] 


V.A 


I 


1 


49 ±11 


1.7 ±3.0 


n c:+3.05 
u -°-0.43 






V.B 


IV 


3 


49 ±7 




i n+0.70 
i - u -0.41 


-0.97 ±0.40 


i n+0.78 
1 - u -0.44 


V.C 


I 


4 


56 ±30 




i n+2.09 
i - u -0.68 




36 ±17 


V.D 


III 


5 


52 ±20 


0.8 ±4.8 


5+ 4 ' 29 


0.01 ±2.62 


0.13±g;f| 35 ±15 



TABLE VII: The 2a predicted measurements for the next generation of ton-scale direct detection experiments 
assuming different DM benchmarks and data sets. In all cases presented here, astrophysical uncertainties are 
marginalised over according to Table |VT| 



only be constrained at the expenses of introducing theoretically motivated hypotheses 
regarding WIMP couplings. 

To conclude, in view of the next generation of ton-scale instruments, it is increasingly 
important to identify what direct detection can really tell us about WIMPs independently 
from any theoretical prejudice (regarding couplings, for instance) and keeping account of 
astrophysical uncertainties. Only that way may one hope to convincingly pinpoint the 
nature of dark matter. 
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